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Resistivity measurements from 4.2 K up to 300 K were made on YNi2B2C single crys- 
tals with T c =15.5 K. The resulting p(T) curve shows a perfect Bloch-Griineisen (BG) 
behavior, with a very small residual resistivity which indicates the low impurity con- 
tent and the high cristallographic quality of the samples. The value Atr =0.53 for the 
transport electron-phonon coupling constant was obtained by using the high-temperature 
constant value of dp/dT and the plasma frequency reported in literature. The BG ex- 
pression for the phononic part of the resistivity p prl (T) was then used to fit the data in 
the whole temperature range, by approximating a^ r F(fl) with the experimental phonon 
spectral density G(Q) multiplied by a two-step weighting function to be determined by 
the fit. The resulting fitting curve perfectly agrees with the experimental points. We 
also solved the real-axis Eliashberg equations in both s- and d-wave symmetries under 
the approximation a 2 F(£l) m a^ r F(^l). We found that the value of Atr here determined 
in single-band approximation is quite compatible with T c and the gap A experimentally 
observed. Finally, we calculated the normalized tunneling conductance, whose compar- 
ison with break-junction tunnel data gives indication of the possible s-wave symmetry 
for the order parameter in YNi2B2C. 



1 Introduction 

In the past seven years, a great interest has been aroused by the discovery of the fam- 
ily of superconducting quaternary rare-earth borocarbide intermetallic compounds 
RNi2B2C (R = rare earth), of which YNi 2 B 2 C is one of the most studied non- 
magnetic members. In some way similarly to high-T c cuprates, the borocarbides 
have a layered crystal structure, even if band-structure calculations have shown 
the presence of a three-dimensional electronic structure (see Ref. [1] for a short 
overview of selected properties). Initially, a conventional BCS description of the 
superconducting properties in these materials has been supported by many of the 
experimental and theoretical results. More recently some particular features of the 
electronic specific heat (as a function of temperature and magnetic field) and of the 
upper critical field in the non-magnetic borocarbides LUM2B2C and YM2B2C have 
been interpreted as signs of an unconventional pairing, possibly of d-wave typeH. 
On the other hand, Shulga et al. have recently explained the positive curvature 
of the upper critical field of YNi2B2C as a function of the temperature near T c , 
its magnitude and shape in the framework of s-wave Migdal-Eliashberg theory by 
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Figure 1: A few points of the typical resistivity of one of our YNi2B2C single crystals (open circles) 
and the fit using the Bloch-Griineisen model (solid line). The inset shows the low temperature 
part of the resistivity and the fit with all the measured points. 

considering the presence of two bands, ona of them being more deeply involved in 
the transport properties of the compound 0. 

In the present work we demonstrate the complete agreement of the experimental 
resistivity data obtained in YM2B2C with the predictions of the classic single- 
band theory for the electron-phonon (e-p) interaction (Bloch-Griineisen theory) . We 
obtain a value of the transport e-p coupling constant that well agrees with previous 
experimental and theoretical results. The transport e-p spectral function a^ T F(fl) 
is also obtained by the fit of the experimental data in the whole temperature range 
and used for the calculation of the normalized tunneling conductance by directly 
solving the real- axis Eliashberg equations. 

2 Experiment 

High-quality YM2B2C single crystals were grown by using the rf - zone melting 
process □. The critical temperature of the crystals, measured by AC susceptibility 
and perfectly confirmed by resistivity measurements, is T c = 15.5 K with AT C (10 — 
90%) = 0.5 K. The imaginary part of the susceptibility shows a single, very high 
and narrow peak at T w T c which confirms the purity and crystallographic quality 
of the samples. 

We measured the resistivity of these crystals as a function of the temperature, 
on cooling from 300 K down to 4.2 K, by using the AC version of the standard four- 
probe technique. The current leads were directly soldered to the opposite sides of 
the samples, which have parallelepipedal shape. The voltage leads, made from very 
thin gold wires, were glued with a conducting paste to the surface of the crystals. 
We improved the sensitivity of the measurement by injecting in the crystals an AC 
current of typically 10 mA at 133 and detecting the voltage by the standard lock- 
in technique. Due to a very slow cooling-down procedure, we were able to collect 
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Figure 2: The relative deviations (p— Pfit)/P of the experimental data from the fit. 

nearly three thousand resistivity values for every curve between 4.2 and 300 K. 

In Figure 1 the resistivity p(T) of one of the YM2B2C crystals is shown (open 
circles). For clarity, only a reduced number of points is reported. In the inset of the 
same figure we show an enlargement of the low-temperature part of the resistivity 
curve that contains all the measured points at T < 30 K (open circles). 

The resistivity of Fig. 1 shows a perfect Bloch-Griineisen (BG) behavior that 
has already been observed in previous experiments El. The small residual value of 
the resistivity p(0)— 3 //il-cm and its perfectly linear high-temperature part (with 
a slope dp/dT = 0.12 pSl-cm/K) indicate the high quality and low impurity content 
of the samples. Quite similar results were obtained in various YM2B2C samples. 

3 Discussion 

The resistivity of a normal Fermi-liquid metal follows the well known Matthiessen's 
rule: p(T)=p +p p i 1 (T), where po and p P h(T) are the residual and the phonon resis- 
tivities, respectively. 

According to the BG theory, the high-temperature part of the p p h(T) shows a 
linear behavior given by the following expression: 

p ph (T) = (2ne k B /hu 2 p )X tI T 

where ui p is the plasma frequency, At r =2 J °° [a% T F(£l) /Sl]dSl is the transport e-p cou- 
pling constant, and a^Fifl) is the transport e-p spectral function. Once known the 
experimental value of the plasma energy hw p and the temperature coefficient of the 
linear part of the resistivity at T>100 K (see Fig. 1), we can determine the transport 
coupling constant Atr = (^- , p/27reofcB)dp p h/dT = (hu>p/2irEokB) dp/dT. The values 
uj p — 4.25 eV determined in Ref. [6] by reflectance and EELS measurements and 
dp/dT=0.12 juQ-cm/K from Fig. 1 lead to the result A tr =0.53. 

Full information on the e-p coupling in YM2B2C can be obtained by fitting the 
resistivity of Fig. 1 in the whole temperature range. Again following the standard 
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Figure 3: The electron-phonon spectral function determined from the fit of the resistivity of Fig. 1. 
In the inset the step-like e-p coupling function obtained from the fit is shown. 

BG theory, we use the most general expression for p p h(T), given by: 

/•n max 

p ph (T) = (4ire k B T/hujl) / [a t 2 r F(ft)/0] • [HSl/2k B T amh(hn/2k B T)] 2 dn. 

Jo 

Our goal is to determine the function a 2 r F(r2) from the fit. Actually, as a first 
approximation, for the a^ T F(fl) we use the phonon spectral density G(f2) obtained 
in inelastic neutron scattering experiments u multiplied by a two-step weighting 
function, whose constant values for 0<37.5 meV and 37.5<f2<70 meV are to be 
determined by the fit. These energy ranges are chosen because they correspond to 
the two most-distinguishable structures of the G(f2) which shows peaks at about 20 
and 50 meV and a value close to zero just at 37.5 meV. We neglected in the fit the 
presence of high energy phonons at about 100 meV. 

The results of the fit are shown as a solid line in Fig. 1. The experimental p(T) 
is perfectly fitted by the theoretical BG curve in the whole temperature range. Fig. 2 
shows the relative deviations of the experimental curve from the fit (p — pfit)/p as 
function of temperature. They never exceed ±5%. 

Fig. 3 shows the spectral function a 2 r F(f2) obtained as the product of the G(fl) 
by the two-step weighting function a 2 r determined from the fit. The latter function 
is shown in the inset of Fig. 3. Of course, the Atr calculated from the af r F(il) 
perfectly coincides with the value previously obtained from the linear part of p{T). 

It is well known that, to a first approximation, Atr ~ A, where A is the e-p 
coupling factor involved in the BCS coupling of the Cooper's pairs. In the hypothesis 
of an e-p coupling origin of the superconductivity in YN12B2C we calculate the 
quasiparticle density of states in this compound in both s- and d-wave symmetry bv. 
solving in direct way the real-axis Eliashberg equations for the strong e-p couplingEl 
and using as a^ r F(fl) the function shown in Fig. 3. 

Fig. 4 shows the SIN tunneling conductances in s- and d-wave symmetry calcu- 
lated at 4.2 K from the density of states (solid and dash lines, respectively). From 
the solution of the Eliashberg equations we obtained both the correct T c and the su- 
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Figure 4: s- and d-wave SIN tunneling conductance of YNi2B2C calculated by direct solution of 
the Eliashberg equations using the e-p spectral function shown in Fig. 3 and the parameters shown 
in the legend. 

perconducting gap A m 2 meV observed in tunneling experimentsi'00, by using a 
value A= 0.57 slightly greater than A tr . This fact is consistent with the conventional 
relation between A and A tr - From the, curves of Fig. 4 and the tunneling experi- 
mental data in STM configurationli3liil it is very difficult to determine the possible 
symmetry of the order parameter in YN12B2C due to the appreciable broadening 
of the SIN data. 

A more clear indication can be obtained by the comparison of the calculated 
SIS tunneling conductance with the break-junction experimental data present in 
literature. This comparison is shown in Fig. 5 where the theoretical SIS curves 
at 4 K determined from the solution of Eliashberg equation for pure d-wave (dash 
line), pure s-jstave (dot line) and s-wave plus magnetic impurities in quasi non- 
unitary limit lla (solid line) are presented together with the break-junction data of 
Ekino et alu (inset). The s-wave curve in presence of a small amount of magnetic 
impurities (for details see the caption of the figure and Ref. [12]) reproduces all the 
main features of the experimental data including the well pronounced dip at about 
twice the energy of the gap. On the other hand, the d-wave tunneling conductance 
is quite different from the experimental curve. These results suggest a possible s- 
wave symmetry (or, at least, a dominant s-wave component) in YM2B2C and give 
evidence for the essential role played by the electron-phonon coupling in the pairing 
mechanism in this compound. 

4 Conclusions 

We have demonstrated that the resistivity of YM2B2C has a temperature depen- 
dence perfectly described by the standard BG theory for the e-p coupling in conven- 
tional metals. The value of At r here determined from the resistivity is representative 
of an intermediate e-p coupling strength and is consistent with the value used in 
Ref. [3] for the discussion of the superconducting properties of YNi 2 B 2 C in the 
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Figure 5: s- and d-wave SIS tunneling conductance of YNi2B2C (dot and dash line, respectively) 
calculated by direct solution of the Eliashberg equations using the e-p spectral function shown 
in Fig. 3 and the parameters indicated in the legend. The solid line represents the s-wave SIS 
conductance for A= 0.57 in presence of a small amount of magnetic impurities corresponding to 
the parameters Co=4 and T=0.8 of the model of Ref. [12]. The inset shows the break-junction 
tunneling conductance at 4.4 K from Ref. [9]. 



framework of the isotropic single-band model. The transport e-p spectral function 
a^F(n) determined from the fit of the resistivity was used, to a first approximation, 
in the direct solution of the Eliashberg equations both in s- and d-w&ve symmetry. 
The resulting T c and A are quite in agreement with the experimental data, while 
the comparison of the calculated SIS tunneling conductance with the break-junction 
data supports a conventional s-wave symmetry for the order parameter in YNi 2 B 2 C. 
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